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PREFACE 


An  analytical  investigation  of  salinity  intrusions  and  related 
phenomena  was  initiated  by  the  Corps  of  Engineers  Committee  on  Tidal  Hy¬ 
draulics  in  January  195^-  The  objectives  of  the  investigation  were  stated 
by  the  Committee  at  its  21st  meeting  at  Vicksburg,  Mississippi,  in 
January  1955: 

"The  present  broad  purpose  of  the  investigation  is  to 
determine  the  effects  of  the  physical  and  hydraulic 
features  of  estuaries  (such  as  tidal  prism,  tidal 
range,  fresh-water  discharge,  channel  depth,  channel 
width,  etc.)  on  the  extent  of  salinity  intrusion,  the 
nature  of  salinity  intrusion,  the  magnitudes  and 
durations  of  current  velocities,  and  other  factors 
considered  essential  to  proper  solution  of  estuarine 
problems  encountered  by  the  Corps  of  Engineers." 

The  experimental  portion  of  the  investigation  was  conducted  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (Civil  Works  Investiga¬ 
tion  Item  CW  8^3)  under  the  immediate  supervision  of  Mr.  Henry  B.  Simmons, 
Chief,  Estuaries  Section.  All  of  the  experimental  data  presented  and 
analyzed  in  this  report  were  obtained  in  the  Waterways  Experiment  Station 
rectangular  salinity  flume  and  were  forwarded  to  the  authors  of  this  re¬ 
port  during  the  period  between  September  1955  and  November  1959*  The 
analytical  part  of  the  study  was  carried  out  by  the  authors.  Dr.  Arthur  T. 
Ippen  and  Dr.  Donald  R.  F.  Harleman,  at  the  Hydrodynamics  Laboratory, 
Massachusetts  Institute  of  Technology,  under  contract  DA-22 -079-eng-l 58. 

The  general  investigation  is  designed  to  cover  the  following  four 
phases:* 


*  Minutes  of  the  3^th  Meeting  of  the  Committee  on  Tidal  Hydraulics, 
June  1959- 
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a.  The  extent  of  salinity  intrusion  and  the  mean  salinity 
distribution. 

b.  The  vertical  mixing  of  fresh  and  salt  water  and  the  result¬ 
ing  vertical  salinity  distribution. 

c.  The  vertical  distribution  of  current  velocities  as  affected 
by  salinity  distribution. 

d.  The  movement  and  deposition  of  sediments  as  affected  by 
density -current  phenomena. 

This  report  is  primarily  concerned  with  the  first  phase,  a  one -dimensional 
treatment  of  salinity  intrusion  in  a  partially  or  well-mixed  type  of 
estuary. 

In  June  1957  the  authors  submitted  a  preliminary  report  to  the  Com¬ 
mittee  on  Tidal  Hydraulics  which  defined  and  verified  the  basic  analytical 
procedures.  This  resulted  in  an  expanded  experimented  program  to  test  the 
conclusions  reached  at  that  time.  The  present  report  should  be  considered 
as  superseding  the  previous  report  in  view  of  the  more  comprehensive  ex¬ 
perimental  program  and  conclusions  now  available. 
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Tidal  amplitude  in  ocean,  ft 

Amplitude  of  incident  and  reflected  waves  at  closed  end  of  estuary,  ft 
Width  of  estuary,  ft 

Seaward  excursion  to  point  where  salinity  is  constant  throughout 
tidal  cycle,  ft 

Tidal  wave  celerity,  ft/sec 

Frictionless  tidal  wave  celerity,  ft/sec 

2 

Turbulent  diffusion  coefficient,  ft  /sec 

2 

Apparent  diffusion  coefficient,  ft  /sec  (subscript  o  indicates 
diffusion  coefficient  at  ocean  entrance) 

Tidal  wave  energy  per  unit  surface  area,  lb/ft 

2 

Acceleration  of  gravity,  ft/sec 

Rate  of  tidal  energy  dissipation  per  unit  mass  of  fluid,  ft  /sec3 
Mean  water  depth,  ft 

2  *3 

Rate  of  gain  of  potential  energy  per  unit  mass  of  fluid,  ft  /sec3 

Wave  number,  2ir/X 

Mean  eddy  size,  ft 

Length  of  estuary,  ft 

Minimum  salinity-intrusion  length,  ft 

Maximum  salinity -intrusion  length,  ft 

Manning’s  roughness 

Energy  flux  (power),  ft-lb/sec 

Tidal  prism,  ft 
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Fresh-water  discharge  into  estuary,  ft ■y sec 
Salinity  concentration,  ppt 
Salinity  concentration  in  ocean,  ppt 
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t  Time,  sec 

t„  Local  time  of  high  water 
n 

t^  Local  time  of  maximum  velocity 
T  Tidal  period,  sec 

u  Tidal  velocity  in  x  direction,  ft/sec 
uq  Maximum  tidal  velocity  at  ocean  entrance,  ft/sec 
U  Velocity  of  fresh -water  flow,  ft/sec 
v  Velocity  in  y  direction,  ft/sec 
w  Velocity  in  z  direction,  ft/sec 
W  U/2D^B 

x  Longitudinal  direction,  ft 
y  Vertical  direction,  ft 
z  Width  direction,  ft 
a  Tan-1  (0/2it) 

P  Dimensionless  time  interval,  t/T 

7  Specific  weight  of  fluid,  lb/ft^ 

AD  Mass-transfer  coefficient  due  to  gravitational  convection 

8  ao/huQ 

T)  Tidal  wave  elevation  above  or  below  mean  sea  level,  ft 
t)h  Local  high  water,  ft 

1  „  High  water  at  closed  end  of  estuary,  ft 
On 

6  at  +  a 
M 

X  Tidal  wave  length,  ft 

Xq  Frictionless  tidal  wave  length,  ft 

p  Frictional  wave  damping  factor,  l/ft 

2/1+ 

p  Fluid  density,  lb-sec  /ft 
a  Tidal  frequency,  2n/T 


ONE-DIMENSIONAL  ANALYSIS  OF  SALINITY  INTRUSION  IN  ESTUARIES 


PART  I:  INTRODUCTION 

1.  The  objective  of  the  investigation  reported  herein  is  to  develop 
analytical  procedures  for  the  determination  of  the  distribution  of  salinity 
in  an  estuary.  The  intrusion  of  saline  water  in  the  tidal  portion  of  a 
river  involves  the  basic  mechanism  of  mass  transfer  by  turbulent  diffusion 
and  mass  transfer  by  convective  currents  which  are  associated  with  tidal 
motion  and  the  presence  of  liquids  of  different  density.  Usually  the  en¬ 
gineer  is  called  upon  to  predict  the  change  in  the  salinity  distribution  in 
a  particular  estuary  which  will  occur  as  a  result  of  a  change  in  estuary 
geometry  or  hydraulic  characteristics.  Common  examples  of  such  changes  are 
the  dredging  of  deeper  navigation  channels  and  changes  in  fresh-water  flow 
rates  due  to  upstream  diversions  or  regulation.  While  such  problems  have 
been  studied  by  engineers  and  oceanographers  for  many  years,  in  almost  all 
instances  the  concern  has  been  with  a  particular  estuary.  In  order  to 
divorce  the  investigation  from  a  particular  geometry  and  to  simplify  the 
problem,  the  present  study  is  confined  to  the  case  of  a  rectangular  estu¬ 
ary  of  constant  mean  depth  and  width. 

Previous  Investigations 

2.  An  evaluation  of  the  state  of  knowledge  of  salinity-intrusion 

phenomena  in  estuaries  was  prepared  in  1950  by  the  Committee  on  Tidal  Hy- 
16* 

draulics.  A  review  of  estuarine  hydrography  has  been  given  by 
13 

Pritchard  in  which  estuaries  are  classified  in  terms  of  the  directions 
of  the  advective  and  nonadvective  transport  processes.  The  estuary  of 
this  study  is  of  the  partially  mixed  class  (type  B  in  reference  13)  in 
which  mixing  occurs  between  the  fresh  water  moving  down  the  estuary  and 
the  saline  water  moving  up  the  estuary. 

3.  Prior  to  1950  the  analysis  of  longitudinal  salinity  distributions 


*  Raised  numerals  refer  to  similarly  numbered  items  in  the  list  of 
references  at  end  of  text. 
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was  based  on  the  concept  of  the  tidal  prism.  The  most  precise  formulation 
of  the  tidal -prism  approach  is  due  to  Ketchum  who  suggested  dividing  an 
estuary  into  segments  having  lengths  equal  to  the  average  distance  tra¬ 
versed  by  a  particle  of  water  on  the  flood  tide.  Within  each  segment  it  is 
assumed  that  there  is  complete  mixing  at  high  tide.  The  method  is,  there¬ 
fore,  only  applicable  for  those  estuaries  which  can  be  said  to  approach  the 
fully  mixed  condition.  Such  estuaries  are  characterized  by  extremely  large 
ratios  of  tidal  prism  to  river  discharge. 

4.  For  example,  even  though  the  Delaware  estuary  has  a  tidal  prism 
to  fresh-water  discharge  ratio  of  the  order  of  100,  it  still  exhibits  char¬ 
acteristics  of  a  partially  mixed  estuary  in  the  sense  that  salinity  trans¬ 
fer  by  gravitational  convection  currents  is  quite  strong. 

5-  Arons  and  Stommel1  have  proposed  a  mixing-length  theory  of  tidal 
mixing  as  a  means  of  treating  the  time  average  (over  a  tidal  cycle)  salin¬ 
ity  distribution  in  a  rectangular  estuary.  The  one -dimensional 
conservation-of-salt  equation  was  employed  with  a  convective  term  for  the 
river  flow  and  a  horizontal  eddy  diffusivity.  The  latter  is  assumed  to  be 
equal  to  the  product  of  the  maximum  tidal  velocity  at  the  estuary  entrance, 
the  tidal  excursion  length,  and  a  constant  of  proportionality.  By  inte¬ 
grating  the  conservation  equations,  a  family  of  salinity  distribution 
curves  is  obtained  in  terms  of  the  distance  along  the  estuary  divided  by 
the  total  length  of  salinity  intrusion.  The  results  are  primarily  useful 
as  a  classification  of  estuaries  by  means  of  a  "flushing  number"  obtained 
by  a  best  fit  of  field  salinity  measurements  with  one  of  the  family  of 
curves . 

6.  By  contrast,  the  present  report  will  treat  the  problem  of  pre¬ 
dicting  the  instantaneous  salinity  distribution,  and  therefore  the  intru¬ 
sion  lengths,  in  terms  of  the  tidal  and  hydraulic  properties  and  certain 
boundary  conditions  to  be  determined  at  the  ocean  entrance  to  the  estuary. 
It  is  therefore  necessary  to  consider  the  time -dependent  conservation-of- 
salt  equation.  It  is  this  approach  which  will  be  followed  in  the  present 
study. 

Experimental  and  Analytical  Program 


7.  The  experimental  work  was  conducted  at  the  Waterways  Experiment 
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Station  in  a  horizontal,  rectangular  lucite  flume  327  ft  long  and  9  in. 
wide.  The  mean  water  depth  in  the  flume  was  6  in.  in  the  majority  of  the 
tests.  One  end  of  the  flume  is  closed  and  contains  provision  for  the  in¬ 
troduction  of  fresh  water  at  a  controlled  rate;  the  opposite  end  of  the 
flume  is  connected  to  a  large  tidal  basin  wherein  a  constant  salinity  and 
a  simple  harmonic  tide  can  be  maintained.  Twenty  runs  representing  various 
test  conditions  were  used  in  the  analysis;  however,  the  total  number  of 
flume  tests  was  larger  than  this  because  of  duplication  of  certain  runs  and 
incomplete  data  in  some  of  the  earlier  tests. 

8.  A  detailed  description  of  the  salinity  flume  and  the  method  of 
conducting  the  tests  is  contained  in  reference  17-  A  summary  of  the 
physical  characteristics  of  the  flume  tests  is  given  in  table  1  in  terms  of 
tidal  amplitude  and  period,  salinity  in  the  tidal  basin,  fresh-water  dis¬ 
charge,  and  flume  roughness.  The  roughness  was  obtained  by  use  of  lA-in.- 
square  lucite  strips,  spaced  on  2-in.  centers.  For  the  roughness  desig¬ 
nated  as  "side,"  the  strips  were  attached  to  the  sidewalls  only;  the  desig¬ 
nation  "bottom"  indicates  that  the  strips  were  attached  to  the  bottom  only; 
and  "smooth"  indicates  the  absence  of  strips. 

9-  In  a  rectangular  estuary  of  constant  cross  section,  with  a  tidal 
range  sufficient  to  promote  vertical  mixing,  the  longitudinal  salinity  dis¬ 
tribution  may  be  treated  one -dimensionally.  In  this  sense  the  term  salin¬ 
ity  s  refers  to  the  average  concentration  of  salt  at  a  particular  sec¬ 
tion  x  and  instant  of  time  t  . 

10.  Using  the  notation  shown  in  fig.  1,  the  one -dimensional  salinity 
distribution  problem  can  be  considered  in  terms  of  nine  independent  vari¬ 
ables.  In  general, 

s  =  /  (x,  t,  L,  h,  a,  T,  U,  sq,  and  roughness) 

TIDAL  PER»OD=  T 


Fig.  1.  Notation  for  one -dimensional  salinity  intrusion 
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11.  The  specific  objective  of  this  phase  of  the  investigation  is  the 
correlation  of  the  experimental  data  by  means  of  an  analytical  study  of  the 
salinity- intrusion  mechanism.  Throughout  the  analysis,  the  ultimate  objec¬ 
tive  of  applying  the  information  gained  in  the  one -dimensional  study  to 
actual  estuaries  is  of  primary  importance.  Whenever  a  choice  of  analyt¬ 
ical  procedures  is  available,  the  one  most  applicable  to  natural  conditions 
is  chosen.  Approximations,  which  in  certain  cases  axe  unnecessary  in  the 
one -dimensional  problem,  are  introduced  (sometimes  at  the  expense  of  mathe¬ 
matical  rigor)  to  further  achieve  this  goal. 
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PART  II:  THEORETICAL  CONSIDERATIONS 

12.  Salinity  intrusion  is  closely  associated  with  the  tidal  motion 
in  an  estuary.  The  tides  are  the  primary  source  of  energy  for  the  turbu¬ 
lent  mixing  of  fresh  and  salt  water;  in  addition,  large  convective  motions 
are  associated  with  the  great  length  of  tidal  waves.  For  convenience,  the 
theoretical  considerations  will  be  presented  in  two  sections,  the  first 
dealing  only  with  the  tidal  dynamics  and  the  second  with  the  intrusion  of 
salt. 


Tidal  Dynamics 

13*  At  the  mouth  of  an  estuary  a  cyclical  rise  and  fall  of  the  water 
surface  exists  as  the  result  of  tidal  action  in  the  ocean.  The  tidal  mo¬ 
tion  within  the  estuary  is  entirely  the  result  of  tide -producing  forces  at 
the  mouth  and  is  termed  a  co-oscillating  tide.  If  friction  is  neglected, 
the  co-oscillating  tide  becomes  a  standing  wave  and  the  extreme  water 
levels  occur  simultaneously  over  the  entire  basin.  A  standing  wave  in  an 
estuary  having  a  length  L  less  than  one-quarter  the  length  of  the  stand¬ 
ing  wave  X  is  shown  in  fig.  2. 


Fig.  2.  Standing  wave  in  a  channel  where  L  <  X/4 


In  this  case  the  range  of  tide  at  the  end  of  the  estuary  is  considerably 
greater  than  that  at  the  ocean  entrance.  For  frictionless  tidal  waves, 

c  =  N/gh  and  T  =  12.4  hr,  X  =  c  T  =  V gh  T  .  For  an  estuary  having  a 

o  00' 

mean  depth  of  50  ft,  the  wave  length  is  325  miles;  thus  the  lengths  of  many 
natural  estuaries  are  less  than  X/4. 
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14.  A  study  of  the  U.  S.  Coast  and  Geodetic  Survey  tide  tables  re¬ 
veals  that  few,  if  any,  natural  estuaries  are  accurately  represented  by  a 
frictionless  standing  wave.  It  is  generally  observed  that  the  time  of 
high  water  in  the  landward  direction  is  progressively  later  with  respect 
to  high  water  at  the  ocean  end.  In  a  uniform  estuary  this  time  lag  must 
be  entirely  attributed  to  frictional  effects.  In  the  damped  co-oscillating 
tide  the  amplitude  of  the  incoming  wave  decreases  in  the  landward  direc¬ 
tion,  is  reflected  at  the  end  of  the  estuary,  and  again  undergoes  damping 
in  its  progress  toward  the  ocean.  In  the  limiting  case  the  incoming  wave 
is  completely  dissipated  by  friction  and  the  amplitude  at  the  head  of  the 
basin  and  hence  that  of  the  reflected  wave  is  zero.  In  this  case  the  co¬ 
oscillating  tide  reduces  to  a  single  progressive  wave  subject  to  extreme 
damping. 

15-  The  basic  theory  of  the  co-oscillating  tide  is  based  on  the 

2  5  4 

classical  work  of  Defant  and  Fjeldstad.  Einstein  and  Fuchs  have  used 

a  similar  approach  to  predict  tidal  amplitudes  and  velocities  in  a  uniform 

rectangular  estuary,  the  hydraulic  roughness  of  which  is  a  known  constant. 

12 

Recently  Perroud  has  extended  this  work  to  include  prismatic  channels 
whose  shape  can  be  described  by  linear  or  exponential  expressions.  These 
investigators  have  linearized  the  equation  of  motion  by  replacing  the 
quadratic  friction  term  by  a  constant  times  the  velocity.  This  is  ra¬ 
tionalized  on  the  basis  that  the  frictional  work  done  over  a  tidal  cycle 
is  made  identical  for  both  the  linear  and  the  quadratic  damping.  However, 
the  problem  of  the  (a  priori)  choice  of  the  proper  numerical  value  of  the 
friction  coefficient  for  the  linear  damping  term  has  not  been  resolved.  A 
necessary  additional  assumption  is  that  this  coefficient  must  be  constant 
with  respect  to  x  .  Stepwise  procedures  for  tidal  computations  have  been 

■5 

outlined  by  Dronkers  and  Schoenfeld;  in  all  cases  bottom  frictional  coef¬ 
ficients  and  tidal  velocities  or  elevations  must  be  known  boundary 
conditions . 

16.  A  somewhat  different  approach  based  on  the  assumption  of  expo¬ 
nential  damping  of  the  incident  and  reflected  waves  was  suggested  by 

14  9 

Redfield  and  was  also  used  by  the  authors  in  the  development  of  a 

mathematical  model  for  the  tidal  dynamics  of  the  Bay  of  Fundy.  The  proce¬ 
dure  requires  that  tidal  ranges  be  known  at  several  points  within  the 
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estuary.  In  the  case  of  a  rectangular  estuary  this  information  permits  the 
calculation  of  the  damping  coefficient  and  the  tidal  velocities.  The  pri¬ 
mary  advantage  of  the  method  is  that  it  can  be  applied  to  natural  estuaries 
when  it  is  desired  to  have  information  on  the  variation  of  tidal  velocities 
along  the  axis  of  the  estuary  for  use  in  salinity-intrusion  studies.  In 
these  cases  the  so-called  "damping  factor"  includes  both  the  effect  of 
cross-sectional  area  changes  and  frictional  resistance.  The  annual  tide 
tables  of  the  Coast  and  Geodetic  Survey  provide  sufficiently  detailed  in¬ 
formation  on  tidal  ranges  in  estuaries  for  the  approximate  calculation  of 
local  velocities.  While  measurements  of  tidal  velocities  in  the  Bay  of 
Fundy  are  scanty,  the  available  data  seem  to  confirm  the  predicted  values. 

17.  Both  the  method  outlined  above  and  the  method  of  Einstein  and 
Fuchs  can  be  applied,  with  equal  facility  and  mathematical  rigor,  to  the 
uniform  rectangular  estuary.  The  Einstein -Fuchs  theory  was  successfully 
compared  with  the  Waterways  Experiment  Station  data  in  the  preliminary  re¬ 
port  by  the  authors  in  June  1957-  In  the  present  report  the  formulation 
described  in  the  preceding  paragraph  will  be  used  in  view  of  its  immediate 
applicability  in  the  case  of  natural  estuaries. 

Co-oscillating  Tide  with  Exponential  Damping 


l8.  Fig.  3  is  a  definition  sketch  showing  the  notation  to  be  used 


in  the  tidal  analysis. 

INC  DENT  WAVE 

MSL  - 

REFLECTED  WAVE 

.  - W' 

OCEAN  X=_L 

*s0  — Jyys 
H 

^  L 

Fig.  3*  Definition  sketch  for  a  damped  co-oscillating  tide  in  a 

uniform  channel 


Tidal  elevations 

19.  The  tide  in  the  channel  is  assumed  to  be  composed  of  two  waves 
traveling  in  opposite  directions:  (a)  an  incident  wave  entering  from  the 
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ocean  (traveling  in  the  positive  x  direction)  and  having  a  surface  eleva¬ 
tion  which  is  a  function  of  x  and  t  (time);  and  (b)  a  reflected 

wave  (traveling  in  the  negative  x  direction)  having  a  surface  elevation 
rig  ,  also  a  function  of  x  and  t  .  If  it  is  assumed  that  the  amount  of 
energy  dissipated  is  proportional  to  the  total  energy  of  the  wave,  the  in¬ 
troduction  of  friction  leads  to  an  exponential  decrease  of  the  amplitude. 
The  surface  elevations  of  the  incident  and  reflected  waves,  respectively, 
are  then  given  by  the  following  equations: 


=  Aq  e"^*  cos  (at  -  kx) 

(1) 

=  A  e^  cos  (at  +  kx) 

0  '  ' 

(2) 

where 

o  =  2n/T  (T  =  period  of  ocean  tide) 
k  =  2n/X  (X  =  wave  length) 

Aq  s  amplitude  of  wave  at  end  of  channel  (x  =  0) 

U  =  damping  coefficient 

t  =  0  (corresponds  to  high  water  at  end  of  channel) 
e  =  2.718 

20.  The  tidal  elevation  tj  (above  or  below  msl)  at  any  time  t  and 
place  x  in  the  channel  is  given  by  the  sum  of  the  amplitudes  of  the  in¬ 
cident  and  reflected  waves;  thus 

q  =  ti^  +  rig  =  Aq  [e-MX  cos  (at  -  kx)  +  e^X  cos  (at  +  kx)  ]  (3) 

The  time  of  high  water  at  any  position  in  the  estuary  occurs  when 
di]/dt  =  0;  therefore 

^  =  e”^X  sin  (at  -  kx)  +  e^*  sin  (at  +  kx)  =  0  (4) 

Equation  4  may  be  expanded  and  regrouped  in  the  following  form: 


(e^*  +  e_|iX)  sin  at  cos  kx  +  (e^X  -  e_|JX)  cos  at  sin  kx  =  0 
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and  since 


e-4X 


+  e 


-px 


tanh  px  , 


tan  ot„  =  -tan  kx  tanh  px 


(5) 


21.  crt  is  the  relative  time  of  high  water  at  any  position  in  the 
£1 

estuary  with  respect  to  the  time  of  high  water  at  the  end  of  the  channel 
(crtjj  =  0).  From  equation  5  the  value  of  ot^  ,  in  degrees,  is  given  by 


ot 


H 


tan 


-1 


( -tan  kx  tanh  px) 


(6) 


In  order  to  determine  the  damping  coefficient  p  from  the  experimental 
tidal  elevation,  it  is  convenient  to  express  the  local  high-water  tidal 
amplitude  t)^  in  terms  of  the  high-water  amplitude  at  the  end  of  the  chan¬ 
nel  .  The  local  high-water  amplitude  at  any  point  along  the  channel 

can  be  obtained  by  substituting  the  time  angle  for  high  water  ( ot  )  from 
equation  6  into  equation  3>  whence 

t)h  =  2Aq^^(cos  2kx  +-  cosh  2px)  (T) 

At  the  end  of  the  estuary  (x  =  0),  cos  (o)  =  cosh  (0)  =1;  therefore, 

=  2Aq  and  the  desired  ratio  of  local  high  water  to  high  water  at  the 
landward  end  is 


=  '\/^{cos  2kx  +■  cosh  2px) 

\H  V  ^ 


(8) 


At  the  ocean  entrance  (x  =  -L),  n  =  a  is  the  amplitude  of  the  ocean 
tide;  hence 


^oH 


cos  2kL  +-  cosh  2pL) 


(9) 


In  a  channel  of  uniform  cross  section  and  roughness,  both  k  and  p 
would  be  expected  to  be  constants.  This  may  be  stated  algebraically  as 


P 


(10) 


and  equations  6  and  8  may  be  rewritten  as  follows: 


at  =  tan 
n 


2a_  _■ 


[- 

f  r 

1  COS 

2  L 


tan  kx  tanh 


M 


cos  2kx  +  cosh 


A  plot  of  the  tidal  amplitude  ratio  rj  /n  versus  the  time  angle  of  high 

ri  oH 

water  at„  may  therefore  be  made  with  the  independent  variables  $  and 
n 

kx  defining  a  family  of  curves.  Thus,  if  the  tidal  amplitude  and  time  of 
high  water  are  known  at  several  points  along  the  channel,  these  values 
when  plotted  in  the  manner  suggested  above  should  fall  along  a  curve  of 
$  =  constant.  In  principle,  both  p  and  k  could  be  determined  by  know¬ 
ing  the  time  and  amplitude  of  the  tide  at  the  closed  end  of  the  channel 
and  at  only  one  other  intermediate  section.  In  practice,  the  graphical 
method  described  above  gives  values  of  p  and  k  which  are  the  best  fit 
for  all  of  the  available  data. 

Tidal  velocities 

22.  The  equation  for  the  tidal  velocity  can  be  derived  by  making 
use  of  the  continuity  equation  for  unsteady,  nonuniform  motion. 

23.  The  continuity  equation  for  the  incident  wave  (considering 
amplitudes  to  be  small  compared  to  depth)  is  given  by: 

^1  iN  „„ 

-  h  "5t  (13) 


where  u^  is  the  velocity  due  to  the  incident  wave  of  amplitude 
With  given  by  equation  1, 

^1 

-=77-  =  -A  e  ^  a  sin  (crt  -  kx) 
ot  o 


and  the  continuity  equation  13  becomes: 


du.  A  a 

-3^  =  ~T—  e_,IX  sin  (at  -  kx) 


u^  is  obtained  by  integrating  with  respect  to  x 
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A  a 

u^  =  —  /  e-|iX  sin  (crt  -  kx)  dx 


and  by  expanding  the  quantity  in  parentheses, 


u^  =  —  (sin  at  /  e-^  cos  kxdx  -  cos  at  /  e-tix  sin  kxdx) 


A  a  ' 

L 


sin  at  e 


f  -u  cos  kx  +  k  sin  kx  A 

V  2,2  ) 

-ux  (  -4  sin  kx  -  k  cos  kx  \ 

-  cos  at  e  l  - 2 - o - )  +  F(t)  (x5) 


-px  (  -4  sin  kx  -  k  cos  kx 

\  2,2 
M  +  k 


The  integration  constant  F(t)  may  be  evaluated  by  introducing  the  boundary 
condition  that  u^  =  0  at  x  =  0  (end  of  basin)  for  all  values  of  time; 
hence,  F(t)  =0  . 

2k.  Equation  15  may  be  considerably  simplified  by  the  following 
operations: 

A  a  e-px 

u  =  - 2 - 2  ^sin  ( -M  cos  ^  +  k  sin  kx) 

4  +  k 

-  cos  at  ( -4  sin  kx  -  k  cos  kx) ] 

A  a  g-px 

u  =  -*rr - 2 - 2  ^Sin  kx  Bin  ^  +  COS  kx  cos  <jt) 

1  \x  +  k 

+  n  (sin  kx  cos  at  -  cos  kx  sin  at)] 


A  a  -4X 

u  =  - g - 2  COS  ^  "  kx)  +  ^  sin  (kx  "  crt)  1 

V  +  K 

_  Ao  g  e-MX  k  cos  (at  -  kx)  +  4  sin  (kx  -  at) 

1  h  \/~2  ~2  \f~2  ~2  \f~2  ~2 

Vu+k  Vn+k  V  4  +  k 


Let  cos  a  = 


and  sin  a  = 


r 2  , 2 

4  +  k 


<2  ~2 
4  +  k 


Then  tan  oc  =  sin  ^  ^  and  therefore  a  =  tan"^  ^  . 

cos  a  k  k 
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From  equation  10,  it  follows  that 


a  =  tan 


\2rr  ) 


Returning  again  to  equation  l6, 

Ao  0  e"1^ 

u  =  — - - -  ■■  [cos  a  cos  (at  -  kx)  +  sin  a  sin  (kx  -  at)] 

1  -\f~2  ~2 

V  4  +  k 

A  a  -qx 

u,  =  - - -  [cos  (kx  -  at  -  a)] 

1  n  yf~2  ~2 

VM  +  k 

However,  since  the  cosine  is  an  even  function, 

A  a  -qx 

u,  =  — - - - -  [cos  ( at  -  kx  +  a)  ]  (l£ 

1  h  ,/-§ — ~2 

Vu  +  k 

25.  By  a  similar  procedure  the  continuity  equation  is  written  for 
the  reflected  wave,  with  the  result: 


Ac  px 

u  =  -  - - — - —  [cos  (at  kx  +  a)]  (19) 

2  h  ,/2 - jT 

Vu  +  k 

26.  The  tidal  velocity  at  any  time  and  position  along  the  channel 
is  given  by  the  sum  of  u^  and  u^  ;  thus, 

A  a  1 

u  _  _J2 - [e-tJJC  cos  (  at  -  kx  +  a)  -  e^*  cos  ( at  +  kx  +  a)](20) 

h  i/l  71 


p +  k 


Since  Aq  =  l^2  and  tjoH  “ay  be  eliminated  by  means  of  equation  9,  the 
above  equation  may  be  rewritten  as 


The  maximum  velocity  at  the  ocean  end  of  the  channel  will  be  of  interest 
in  the  salinity  analysis;  the  maximum  velocity  at  any  position  occurs 
when  du/dt  =  0;  therefore,  from  equation  21: 
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otj^  =  tan”^"  (-tan  kx  coth  px)  -  cl 


This  value  of  ot^  ,  when  substituted  into  equation  21,  will  give  the 
maximum  velocity  at  any  position  x  .  By  putting  x  =  -L  in  equations 
21  and  22,  the  maximum  velocity  at  the  ocean  entrance  can  be  calculated. 
Tidal  energy  dissipation 

27-  One  of  the  important  concepts  in  turbulent-diffusion  theory  is 
the  relation  between  the  diffusion  coefficient,  the  rate  of  energy  dissipa¬ 
tion  per  unit  mass  of  fluid,  and  the  scale  of  turbulence.  In  a  partially 
or  well -mixed  estuary,  the  primary  source  of  mixing  energy  is  the  tidal  mo¬ 
tion.  The  rate  of  energy  dissipation  per  unit  mass  due  to  the  tidal  motion 
should  therefore  be  a  useful  factor  in  the  correlation  of  the  salinity  in¬ 
trusion  with  the  tidal  characteristics  of  the  channel. 

28.  For  waves  of  small  steepness,  the  energy  of  a  progressive  os¬ 
cillatory  wave  per  unit  of  surface  area  is 


where  ax  is  the  local  maximum  amplitude  of  the  wave.  The  flux  of  en¬ 
ergy  P  across  a  vertical  section  is  given  by 


P  =  Ecb  (ft~—) 
'  sec 


where  c  is  the  wave  celerity  and  b  the  channel  width. 

29.  In  the  damped  co-oscillating  tidal  analysis,  the  tidal  motion 
is  determined  by  the  superposition  of  two  waves,  one  (the  incident  wave) 
proceeding  from  the  ocean  entrance  to  the  closed  end  of  the  estuary  and 
the  other  (the  reflected  wave)  from  the  closed  end  to  the  entrance.  Thus, 
the  energy  flux  for  the  incident  wave  (moving  in  the  +x  direction)  of 
amplitude  kQ  e  ***  (see  equation  l)  is 


,  _  A2  p-M* 

1  2  o  ' 


and  that  for  the  reflected  wave  (moving  in  the  -x  direction)  is 


w 
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=  _  cb21  a2  ^2 px 


2  Ao  6 


The  total  energy  flux  at  an  intermediate  section  is  therefore 


£|*  A2  (e24X  -  e'2^) 


or 

P  =  -cb?  A2  sinh  2px  (23) 

(Here  it  is  seen  that  if  p  =  0  (standing  wave)  the  total  energy  flux 

P  =  0.) 

30.  The  energy  flux  at  the  ocean  entrance  (x  =  -L)  is  therefore 

2 

P  =  cby  A  sinh  2pL 
L  O 


The  change  in  the  energy  flux  between  the  ocean  entrance  and  any  inter¬ 
mediate  section  is  equal  to  the  rate  of  energy  dissipation  in  that  portion 
of  the  channel ;  hence. 


Rate  of  energy  dissipation  =  P 

L 


P  =  Cbr  A 


2 

o 


(sinh  2pL  +  sinh  2px) 


The  mass  of  fluid  in  this  same  portion  of  the  channel  is  pbh  (L  +  x) 
and  the  rate  of  energy  dissipation  per  unit  mass  Gx  is  given  by 


G  ££.  ( sinh  2pL  +  sinh  2px 

x  "  h  o  \  L  +  x 


(24) 


The  rate  of  energy  dissipation  per  unit  mass  for  the  entire  channel  is 
given  by  setting  x  =  0  ;  hence. 


sinh  2pL 


(25) 


„  ,  sinh  2ux 

G  1  +  .  r  o  T 

x  sinh  2pL 


The  ratio 
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and  the  maximum  value  of  this  ratio  occurs  at  the  ocean  entrance  where 


x  — ►  -L  ,  in  the  limit  the  ratio  r  •  If  the  argument  of  the 

’  G  tanh  2pL 


hyperbolic  function  (2pL)  is  small,  the  maximum  value  of  Gx/G  will  not 
differ  appreciably  from  unity;  hence,  it  is  generally  sufficient  to  use 
the  total  energy  dissipation  rate  as  given  by  equation  25  as  representa¬ 
tive  of  the  entire  channel.  The  factor  A  can  be  eliminated  by  means  of 

o 

equation  9,  and  equation  25  thus  becomes 


G  = 


gc_a 

2hL 


sinh  2uL 


cos  2kL  +  cosh 


(26) 


Note  that  in  the  above  equation  the  wave  celerity  c  =  X/T  =  a/k  .  In 
general,  due  to  the  effect  of  friction,  the  celerity  will  be  less  than 
■</gh  as  in  the  frictionless  case. 

Summary 

31.  For  convenient  reference  in  the  later  sections,  the  important 
equations  for  the  tidal  dynamics  are  summarized  below. 

32.  The  two  independent  parameters  describing  the  tidal  motion  are 
the  wave  length  expressed  by  k  =  2?r/X  and  the  damping  modulus  p  .  If 
the  amplitude  and  time  of  high  water  at  the  closed  end  and  at  various  sta¬ 
tions  along  the  channel  are  known,  p  and  k  can  be  determined  from  the 
following  two  equations: 


ot  =  tan 
n 


-i 


(-tan  kx  tanh  px) 


<rtH 


=  0  at  closed  end) 


cos  2kx  +  cosh 


2px) 


(6) 


(8) 


( T)  =  elevation  of  high  water 
at  closed  end) 


The  maximum  tidal  velocity  at  the  ocean  entrance  of  the  channel  uq 
given  by  equation  21  with  x  =  -L  , 


is 


‘i 
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where 


u 

o 


a g  cos  ( S  +  kL)  -  e~^L  cos  (6  -  kL) 


kh 


"^2 ^ j  +  1  (cos  2kL  +  cosh  2pL) 


(21a) 


0  = 


+ 


a  =  tan 


-1 


( -tan  kL  coth  pL) 


(22a) 


The  rate  of  energy  dissipation  per  unit  mass  of  fluid  due  to  tidal  fric¬ 
tion  is  given  by: 


G 


_  gg  a 


c 


sinh  2pL 


2khL  \cos  2kL  +  cosh  2pL 


) 


(26) 


Salinity  Intrusion 


33*  The  basic  relation  in  the  salinity-intrusion  analysis  is  the 
equation  for  the  conservation  of  the  mass  of  salt,  otherwise  known  as  the 
convective-diffusion  equation.  For  a  turbulent,  incompressible  fluid  this 
equation  has  the  following  general  form  (see  reference  15,  page  159): 


where  u  ,  v  ,  w  axe  the  x  ,  y  ,  z  components  of  the  fluid  velocity  and 

2 

D  ,  D  ,  D  are  turbulent  diffusion  coefficients  (ft  /sec), 
x  y  z  N  ' 

3^.  The  above  equation  expresses  the  condition  that  the  net  mass 
transfer  of  salt  by  convection  and  diffusion  into  an  elemental  volume  is 
equal  to  the  time  rate  of  accumulation  of  material  in  that  volume. 

35*  In  the  one -dimensional  analysis,  the  salinity  gradients  in  the 
y  and  z  directions  are  neglected  and  equation  27  becomes 

36.  At  any  point  in  the  channel,  the  x  component  of  the  fluid  ve¬ 
locity  is  equal  to  the  sum  of  the  velocity  due  to  tidal  motion  u(x, t) 
and  the  fresh-water  velocity  -U  j  hence, 


IT 


n  ♦«**)  I  - « I  - 1  (»„  %i)  <®> 

37*  In  the  salinity-intrusion  analysis  it  will  be  convenient  to 
take  x  =  0  at  the  ocean  entrance  of  the  estuary  (see  fig.  l);  in  addi¬ 
tion,  time  t  =  0  will  refer  to  low  tide  in  the  ocean.  The  difficulties 
associated  with  finding  a  general  solution  for  equation  29  are  overcome  by 
treating  the  problem  in  two  distinct  parts. 

Salinity  distribution  at 
low  tide  (quasi-steady  state) 

38.  It  is  recognized  that  although  the  salinity  concentration  de¬ 
pends  on  time,  it  varies  in  a  periodic  manner  due  to  uhe  tide.  Equation 
29  may  therefore  be  considerably  simplified  by  considering  only  average 
values  where  the  averaging  is  over  a  period  of  time  of  the  order  of  a 
tidal  cycle.  Since  the  tidal  velocity  u(x, t)  is  harmonic,  its  average 
value  over  a  cycle  is  zero  and  the  quasi-steady-state  conservation  equa¬ 
tion  may  be  stated  as 

-u  si =  It  (dx  §f)  (30) 

39-  The  basic  assumption  that  the  quasi-steady-state  salinity  dis¬ 
tribution  can  be  represented  by  the  one -dimensional  convective-diffusion 
equation  30  warrants  further  consideration.  This  equation  represents  an 
equilibrium  between  the  seaward  mass  transfer  by  the  average  fresh-water 
velocity  U  and  the  landward  mass  transfer  due  to  turbulent  diffusion. 
However,  in  the  estuary  intrusion  problem  there  is  an  additional  mechanism 
for  the  transport  of  salt  in  the  landward  direction.  This  process,  which 
will  be  called  gravitational  convection,  is  caused  by  the  difference  in 
specific  weight  of  the  ocean  and  fresh  water.  The  saline  wedge,  which  oc¬ 
curs  in  estuaries  entering  an  ocean  with  a  negligible  tidal  range,  is  an 
extreme  example  of  the  upstream  penetration  of  salt  by  gravitational  con¬ 
vection.  In  the  treatment  of  the  wedge  problem  it  is  possible  to  neglect 
the  effect  of  the  turbulent  diffusion.  However,  even  in  the  so-called 
well-mixed  estuary,  it  is  by  no  means  obvious  that  the  effect  of  gravita¬ 
tional  convection  may  be  neglected.  The  gravitational  convection  of  salt 
depends  on  the  vertical  density  gradients  which  cannot  be  represented  in 
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the  one-dimensional  equation  30*  The  only  possible  basis  for  a  continua¬ 
tion  of  the  one-dimensional  approach  in  the  light  of  this  difficulty  is 
the  assumption  that  the  diffusion  coefficient  can  reflect  both  eddy  dif¬ 
fusion  and  gravitational  convection.  This  assumption  has  been  verified 
by  a  series  of  experiments  conducted  under  the  authors '  supervision  at  the 
Massachusetts  Institute  of  Technology  Hydrodynamics  Laboratory  (see  refer¬ 
ences  6  and  10) . 

40.  The  MIT  tests  were  carried  out  in  a  rectangular  channel  in 
which  homogeneous  turbulence  was  generated  by  the  vertical  oscillation  of 
metal  screens  immersed  in  the  fluid.  The  test  program  included  unsteady - 
and  steady-state  diffusion  and  convection  of  fluids  of  different  densities 
for  various  intensities  of  fluid  turbulence.  It  was  found  that  longitudi¬ 
nal  salinity  distributions  could  be  described  by  the  one -dimensional 
conservation-of-mass  equation  in  which  the  turbulent  diffusion  coefficient 
was  replaced  by  a  larger  apparent  diffusion  coefficient  .  In  other 
words, 

Df  =  D  +  AD  (31) 

XXX 

where  D  is  the  eddy  diffusivity  and  AD  is  a  mass-transfer  coefficient 

X  X 

due  to  gravitational  convection.  These  observations  are  not  unreasonable 
if  it  is  considered  that  the  gravitational  convection  gives  rise  to  large- 
scale  internal  circulations  which  may  be  thought  of  as  eddies.  In  this 
sense  AD^  is  also  a  turbulent  diffusion  coefficient  due  to  eddies  having 
mixing  lengths  many  orders  of  magnitude  larger  than  the  normal  eddy  scale 
of  the  shear -generated  turbulence.  The  apparent  diffusion  coefficient 
D'  should  therefore  replace  the  turbulent  coefficient  D  in  the  equa- 

X  X 

tions  developed  in  paragraphs  33-36-  Equation  30  may  now  be  integrated 
with  respect  to  x  ,  noting  that  the  fresh-water  velocity  U  is  constant 
in  a  uniform  rectangular  estuary. 

=  "Dx  +  F(t)  (32) 

where  F(t)  is  the  integration  constant  for  the  partial  differential 
equation.  F(t)  is  arbitrarily  set  equal  to  zero  since  it  can  be  shown 
that  any  other  value  has  the  effect  of  introducing  a  constant  displacement 


of  the  salinity  in  the  x  direction.  In  so  doing  it  is  assumed  that  the 
salinity  distribution  given  by  equation  32  represents  conditions  at  suc¬ 
cessive  times  of  low  tide  in  the  estuary,  since  at  this  time  the  salinity 
extends  a  minimum  distance  into  the  estuary;  hence, 


Us  =  -D*  I2- 
x  dx 


(33) 


Separating  the  variables  and  integrating  again, 

/  =  -u  /  It  or  £nB  +  cx  =  -u  /  5T  (34) 

s  x  x 

In  order  to  complete  equation  34  it  is  necessary  to  specify  the  manner  in 
which  varies  with  x  .  Further  consideration  of  this  p^int  will  be 

given  following  the  discussion  of  the  second  part  of  the  solution  of  the 
salt-transfer  equation  29. 

Salinity  convection 
within  a  tidal  cycle 

41.  If  it  is  accepted  that  the  salinity  distribution  is  known  at 
low  tide  (from  the  above  analysis),  the  remainder  of  the  one-dimensional 
salinity-intrusion  problem  consists  of  determining  the  salinity  distribu¬ 
tion  at  other  instants  of  time  within  the  tidal  cycle. 

42.  A  physical  appreciation  of  the  reasoning  involved  may  be  gained 
by  imagining  the  channel  to  have  a  vertical  barrier  separating  regions  of 
fresh  and  salt  water  just  prior  to  the  start  of  a  test.  As  the  barrier  is 
removed  the  tide  generator  at  the  ocean  end  is  started  and  fresh  water  is 
introduced  at  the  closed  eri  of  the  estuary.  During  the  first  few  tidal 
cycles  a  large  amount  of  diffusion  between  fresh  and  salt  water  must  take 
place  within  each  cycle.  It  is  not  difficult  to  imagine  that  the  amount 
of  diffusion  taking  place  during  each  tidal  cycle  will  progressively 
diminish  and  that  after  a  large  number  of  cycles  (under  constant  tidal  and 
fresh-water  inflow  conditions)  the  diffusion  during  one  cycle  will  be 
small.  It  is  therefore  postulated  that  an  equilibrium  condition  or  a 
quasi-steady  state  has  been  reached  in  which  the  salinity  distribution 
does  not  change  appreciably  in  a  tidal  cycle  due  to  longitudinal  diffusion. 

The  principal  change  in  salinity  at  a  fixed  point  is  therefore  a  result  of 

i 

i 


the  tidal  velocity  carrying  the  salinity  up  and  down  the  estuary.  The 
diffusion  process  due  to  fresh-water  inflow  and  the  tidal  velocities  sure 
very  important  in  determining  the  form  of  the  salinity  distribution  but 
during  one  tidal  cycle  their  effects  are  considered  negligible.  Qualita¬ 
tively,  this  is  equivalent  to  setting  the  right-hand  side  of  equation  29 
equal  to  zero.  Therefore,  if  the  low-tide  salinity  distribution  is  known, 
the  salinity  distribution  at  any  other  time  in  the  tidal  cycle  may  be 
found  by  solving  the  convection  equation: 

it  +  W1'*)  - u]  It =  0  (35) 

In  partially  and  well -mixed  estuaries  the  tidal  velocities  u(x, t)  are 
generally  large  compared  to  the  fresh-water  velocity  U  ;  hence  equation 
35  may  be  represented  by 


The  general  solution  of  this  form  of  partial  differential  equation  (refer¬ 
ence  8,  page  37l)  is  of  the  form 


x2  =  f  (xx) 

where  X^(s,x, t)  =  and  X2(s,x, t)  =  C2  are  any  two  independent  solu¬ 
tions  of 


dt  = 


dx 


u(x,t) 


ds 

0 


(37) 


One  solution  is  ds  =  0  and  hence  =  s  ;  the  second  solution  can  only 
be  obtained  readily  if  the  variables  in  the  expression  for  the  tidal  ve¬ 
locity  u(x, t)  are  separable.  At  this  point  it  would  be  desirable  to  in¬ 
troduce  the  general  expression  for  the  tidal  velocity,  equation  21,  previ¬ 
ously  developed.  However,  the  x  and  t  variables  are  clearly  not 
separable  in  this  equation,  and  it  is  necessary  to  find  a  suitable  ap¬ 
proximation  for  u(x, t)  which  is  valid  at  least  over  the  region  of 
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salinity  intrusion  in  the  channel.  For  this  purpose  the  continuity  equa¬ 
tion  13  is  written  for  the  tidal  velocity. 


cki  1  dn 

S  -  -  h  at 


(38) 


If  it  is  assumed  that  the  water  surface  in  the  estuary  oscillates  hori¬ 
zontally,  t)  =  -a  cos  at  (t  =  0  at  low  tide).  Then  gj  *  -  ^  sin  at  . 

Upon  integrating,  u(x, t)  =  -  — jp  sin  at  +  F(t)  .  Since  in  this  simple 

case  the  velocities  are  90  degrees  out  of  phase  with  the  elevations,  the 
following  boundary  condition  may  be  used  to  determine  F(t):  u  =  uq  sin  crt 
at  x  =  0  (ocean  entrance).  Therefore 


u(x, t)  =  uq  (l  -  6x)  sin  at 


(39) 


where  5  =  r — 
hu 

o 

This  development  leads  to  a  linear  decrease  of  the  tidal  velocity  in  the 
landward  direction  in  which  the  x  and  t  variables  in  the  equation  are 
separable.  The  point  to  stress  is  that  in  this  process  the  maximum  veloc¬ 
ity  at  the  ocean  entrance  of  the  channel  is  not  specified  and  the  correct 
value  computed  on  the  basis  of  the  damped  co -oscillating  tide  may  there¬ 
fore  be  used  (equation  21a).  A  numerical  comparison  of  equations  39  and 
21a  given  in  a  later  section  shows  that  the  linear  approximation  does  not 
lead  to  serious  errors.  Thus  it  is  possible  to  approximate  surface  eleva¬ 
tions  and  velocity  variations  by  linear  laws  for  damped  tides. 

43.  The  second  solution  of  equation  37 

..  _  dx 
"  u(x,t) 

is  readily  obtained  by  substitution  of  equation  39;  therefore, 

u0  /  sin  «  at  .  /  ^ 

C2  .  (1  -  8*)  e-(a/h)  009  011 


and 
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Therefore  the  general  solution  of  the  salinity-convection  equation  36 
is  given  by  the  following: 


or 


C1  -  f  <C2> 


f-' 

o 


(1  -  8x)  e-(a/h)  COS 


(40) 


If  the  salinity  distribution  is  known  at  t  =  0  ,  low  tide,  the  salinity 
distribution  at  any  other  instant  of  time  can  be  determined  from 
equation  39 • 


Variation  of  the  diffusion 
coefficient  along  the  channel 

44.  A  solution  of  the  quasi-steady-state  salinity  distribution, 
equation  34,  depends  upon  the  assumed  variation  of  the  apparent  diffusion 
coefficient  with  x  . 

a.  The  simplest  assumption  is  that  =  constant  =  D'  . 

Hence,  equation  34  becomes 

-0ns  +  =  -  gr  (41) 


The  boundary  condition  for  the  evaluation  of  must  ex¬ 

press  the  fact  that  the  salinity  is  known  at  some  x  ; 
however,  the  only  fixed  salinity  is  that  in  the  ocean.  It 
is  physically  obvious  that  during  the  low-tide  cycle  the 
salinity  at  the  ocean- entrance  of  the  channel  (x  =  0)  can¬ 
not  be  that  of  the  ocean  itself  due  to  the  dilution  effect 
of  the  fresh-water  flow.  To  circumvent  this  difficulty  a 
model  is  postulated  in  which  the  channel  is  extended  sea¬ 
ward  an  imaginary  distance  B  to  a  point  where  the  ocean 
salinity  s  is  maintained  constant  at  all  times.  Hence, 
s  =  sq  at  x  =  -B  ,  and  equation  4l  becomes 


s 

s 

o 


-  (x  +  B) 


(42) 


which  expresses  the  low-tide  salinity  distribution  as  a 
simple  exponential  curve. 

b.  A  second  and  more  realistic  assumption  regarding  the  ap¬ 
parent  diffusion  coefficient  (D^  =  E>x  +  ADx)  is  that  it 
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decreases  in  the  landward  direction.  The  reasons  for  the 
decrease  are  twofold: 

(1)  The  eddy  diffusivity  Dx  would  be  expected  to  decrease 
in  the  landward  direction  due  to  the  decrease  in  the 
magnitude  of  the  tidal  velocities. 

(2)  Experiments  in  the  MIT  diffusion  flume"*"  have  shown  that 
ad  is  proportional  to  the  local  density  difference 
along  the  estuary.  Since  Ap/p  decreases  as  x  in¬ 
creases  (in  the  landward  direction),  AD  is  in¬ 
versely  proportional  to  x  . 

The  simplest  assumption  compatible  with  the  above  is  that 
is  an  inverse  function  of  x  having  the  form 


45. 


D* 


x 


C 

x  +  B 


D'  B 
o 

x  +  B 


(43) 


Therefore  at  x  =  0 

D'  =  D' 
x  o 

while  at  x  =  -B 

D'  =  00 

x 

and  for  large  x 

D'  — *~0 
x 

The  fact  that  Dx  approaches  an  infinite  value  at  x  =  -B 
is  consistent  with  the  fact  that  if  the  walls  of  the  chan¬ 
nel  are  imagined  to  extend  to  x  =  -B  an  infinite  amount 
of  mixing  would  be  required  to  maintain  constant  salinity 
throughout  the  tidal  cycle. 


Equation  34  can  be  written  as 


2 

+■  C2  =  -  /  (x  +  B)  dx  -  -  (44) 

o  o 

and  using  the  same  boundary  condition  s  =  sq  at  x  =  -B,  equation  44 
becomes 

_s _ -W(x  +  B)2 

s 

o 


(*5) 


w  = 


u 

2D'  B 
o 


where 
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Equation  45  has  the  form  of  a  Gaussian  or  normal  error  curve.  Fig.  4 
shows  schematically  the  boundary  condition  and  the  characteristic  shape  of 
the  low -water  salinity  distribution  curves  for  the  two  cases  derived  above. 


ELEVATION 

Fig.  4.  Schematic  salinity  distribution  at  low  tide 


The  experimental  results  presented  later  in  this  report  indicate  that  the 
assumption  of  as  an  inverse  function  of  x  is  valid,  and  the  re¬ 

mainder  of  the  analytical  development  is  therefore  based  on  this  assumption. 
Instantaneous  salinity  distributions 

46.  The  preceding  developments  may  be  summarized  as  follows: 

Equation  45  describes  the  salinity  distribution  at  low  tide  (t  =  0) 

-5-  =  e'W(x4  +  B)2  for  (x.  +  B)  >  0  (45) 

s  ii 

o 

As  a  reminder  that  when  x  is  varied  in  equation  45  it  produces  the 


salinity  distribution  at  low  tide  only,  the  subscript  L  is  added  to  the 
x  in  this  equation. 

k'J.  The  salinity  distribution  at  any  other  instant  of  time  is  given 
by  equation  40 


(1  -  Sx)  e"(a/h)  COS  at 


m 


These  equations  (40  and  45)  may  be  combined  to  produce  an  analytical  ex¬ 
pression  for  the  instantaneous  salinity  distribution. 

48.  Assume  that  it  is  desired  to  find  the  distribution  of  salinity 
at  time  t  ;  equation  40  implies  that  salinity  at  a  particular  x  and  t 
will  be  equal  to  the  salinity  at  low  tide  (x^  and  t  =  0) .  The  relation 

between  x  ,  x  ,  and  t  can  be  found  as  follows: 

Lt 

f  [(equation  40)]  =  f  [(equation  40)] 

XT  x 


hence, 

(1  -  8x  )  e-(a/h)  -  (1  -  Sx)  e-(a/h)  cos  at 
h 


Solving  for  x^  : 


=  1  _  (1  -  Sx)  e(a/h)(l  -  cos  0t) 


(46) 


L  8 


Upon  substitution  of  x  into  equation  45  the  instantaneous  salinity  dis- 

Li 

tribution  equation  is  obtained: 


Si***!  =  exp  {.  -V  [i  .  (1  -  8x)  e(a/h)(l  -  COS  rt)+  8b12} 

So  L  2D’  B  8  L  J  J 

o 

The  specification  that  xT  +  B  >  0  in  equation  45  requires  that 

If 

[l  .  (1  .  8x)  e<*/h><l  '  “s  ,  sb]  >  0 


(47) 


in  equation  kj.  If  the  term  in  brackets  is  less  than  zero,  the  ratio 
s(x, t)/so  is  to  be  taken  as  unity. 
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49*  The  maximum  intrusion  of  salinity  occurs  at  high  tide  when 
2jffc 

crt  =  =  n  .  For  this  condition,  equation  47  becomes 


exp  /-—a— 5  [I  .  (1  .  &x)  e2^  +  SB]2}  (W) 

So  l  2D'  B  8^  J 

o 

Again,  sjj^/s0  =  1  if  [l  -  (l  -  8x)  e28^*1  +  8B]  <0  .  The  maximum  or 

total  intrusion  length  L  is  arbitrarily  defined  as  the  point  on  the 

high-tide  salinity  curve  at  which  the  salinity  is  l/lOO  of  the  ocean 

/  -k.6 

salinity-  Hence,  x  =  L ,  when  sTni/s  =  0.01  ,  and  since  e  '  =  0.01, 

7  t  HW  o 


- - - p  [1  -  (1  -  BL. )  e28-/*1  +  SB]2  *  4.6 

2D'  B  8  Z 

o 

Solving  for  L.  with  8  =  ao/hu 
t  o 


The  first  term  in  equation  49  represents  the  convective  intrusion  of  salin¬ 
ity  due  to  the  tidal  excursion,  while  the  second  term  represents  that  part 
of  the  salinity  intrusion  due  to  the  diffusion  and  gravitational  convec¬ 
tion  processes. 

50.  Following  a  similar  approach,  the  minimum  intrusion  length 

is  obtained  directly  from  equation  45  for  the  condition  —  =  0.01  . 

so 


L 


m 


(50) 


Hence,  the  maximum  intrusion  length  can  also  be  written  as 

L+  =  —2-  (l  -  e-2a/h)  +  e"2a/h  L 
t  acr  '  '  m 


(51) 


The  instantaneous  salinity  distributions  and  the  intrusion  lengths  can 
therefore  be  determined  if  the  quantities  D^  and  B  are  known.  The 
following  sections  show  how  they  are  determined  from  the  experimental  data 
and  how  they  are  correlated  with  the  independent  quantities  governing  the 
tidal  and  fresh-water  motions. 


2? 


PART  III:  ANALYSIS  OF  EXPERIMENTAL  RESULTS 

51.  The  analysis  of  the  experimental  results  is  presented  in  two 

sections.  The  first  deals  with  the  determination  of  the  tidal  parameters 
p  (damping  modulus)  and  k  =  2n/wave  length,  and  the  determination  of  the 
maximum  entrance  velocity  uq  .  The  second  section  treats  the  determina¬ 
tion  of  the  apparent  diffusion  coefficient  and  the  boundary  condition 

distance  B  .  A  general  siimmary  of  the  physical  characteristics  of  all  of 
the  experimental  tests  is  presented  in  table  1. 

Determination  of  Tidal  Parameters  p  ,  k  ,  and  u 
_  *  ’ _ o 

52.  Tests  28,  29,  and  30  were  run  specifically  to  obtain  tidal  data 
on  surface  elevations  and  velocities  in  the  WES  flume  without  the  effect  of 
salinity  variations  or  fresh-water  velocity.  The  tidal  parameters  p  and 
k  determined  from  these  tests  may  then  be  compared  with  similar  values  ob¬ 
tained  for  other  tests  to  determine  the  effect  of  varying  amounts  of  fresh¬ 
water  discharge,  ocean  salinity,  and 
flume  roughness. 

53*  The  objectives  of  the  tidal 
analysis  are: 

a.  To  determine  p  and 
k  from  the  experimental 
data  on  water -surface 
elevations . 

b.  To  use  these  values  to 
calculate  the  maximum  ve¬ 
locity  u0  at  the  flume 
entrance  for  use  in  the 
salinity  convection  and 
intrusion  equations. 

£.  To  determine  the  rate  of 
energy  dissipation  in  the 
channel  for  use  as  a  cor¬ 
relating  parameter. 

54.  As  an  example  of  the  calcula¬ 
tion  procedure,  the  analysis  for  test  29 
will  be  presented  in  detail.  Fig.  5  is 
a  plot  of  surface  elevation  rj  (above 
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msl)  versus  time  for  each  station  recorded.  By  graphical  interpolation 
the  amplitudes  and  times  of  both  high  water  and  low  water  are  recorded  for 
each  station,  see  columns  (2)  and  (5)  in  table  2.  The  averages  of  the 
elevations  of  both  high  and  low  water  are  listed  in  column  (3),  since  the 
theory  assumes  a  symmetrical  tidal  wave.  The  ratio  of  the  average  local 
amplitude  to  the  amplitude  at  the  closed  end  r^/ is  then  computed  as 
shown  in  column  (4).  A  similar  averaging  of  the  time  of  high  and  low 
water  is  given  in  column  (6).  Finally  the  time  of  high  water  referred  to 
zero  at  the  closed  end  is  shown  in  columns  (7)  and  (8).  The  values  tabu¬ 
lated  in  columns  (4)  and  (8)  are  plotted  as  shown  in  fig.  6,  and  the  value 


Fig.  6.  Determination  of  p  and  k  from  tidal  amplitudes 
and  time  of  high  water 

of  0  (in  equations  11  and  12)  which  best  fits  the  experimental  data  is 
found  by  trial,  (in  test  2 9,  0  =  2.75 •)  The  value  of  kx  for  each  ex¬ 
perimental  point  along  the  0  =  constant  curve  is  interpolated  from  fig.  6 
as  shown  in  column  (9)  of  table  2.  Since  the  x  distance  (column  10)  is 
known  for  each  point,  the  corresponding  value  of  k  is  found  as  shown  in 
column  (ll)  and  the  average  value  is  used  for  the  test.  In  test  29  this 
average  value  is  k  =  0.186  degree/ft  =  0.0033  radian/ ft.  From  equation  10, 
p  =  ^  k  =  (O.OO33)  =  0.0014  radian/ft.  The  maximum  value  of  the 

tidal  velocity  at  the  ocean  entrance  as  calculated  from  equations  21a  and 
22a  is  u0  =  0.44  ft/sec.  Maximum  velocities  at  other  stations  along  the 
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channel  were  also  calculated  In  the  same  manner.  The  variation  of  u 

max 

along  the  channel,  as  determined  from  the  damped  co-oscillating  tidal 
theory,  is  shown  as  the  solid  curve 
in  fig.  7-  Also  shown  is  the  linear 
velocity  approximation  (equation 
39)j  within  the  region  of  salinity 
intrusion  the  two  curves  agree 
quite  well.  For  test  29  vertical 
velocity  distributions  were  meas¬ 
ured  on  the  center  line  of  the  flume 
at  stations  80,  and  160  ft  from 
the  ocean  entrance.  A  comparison  of 
the  experimental  velocities  at  these 
stations  is  shown  in  fig.  7-  The 
points  plotted  are  88  per  cent  of  the  maximum  average  vertical  center-line 
velocities.  This  factor  accounts  for  the  velocity  variation  in  the  trans¬ 
verse  direction  and  was  determined  by  averaging  velocity  measurements  in 
some  of  the  preliminary  flume  tests. 

55.  Fig.  6  also  shows  the  data  from  several  other  runs  in  series  I 
(tidal  amplitude  =  0.05  ft)  and  series  III  (tidal  amplitude  =  0.10  ft). 

All  of  the  pertinent  tidal  data  are  summarized  in  table  3»  For  all  prac¬ 
tical  purposes  p  and  k  may  be  considered  to  be  quite  independent  of 
changes  in  salinity  and  fresh-water  velocity.  The  damping  factor  p  in¬ 
creases  with  increasing  tidal  amplitude  and  decreasing  tidal  period  (test 
30-3) •  These  changes  in  p  due  to  changes  in  tidal  amplitude  and  period 
are  attributed  to  the  assumption  of  linear  frictional  damping. 

56.  The  roughness  indicated  in  table  3  as  side,  bottom,  and  smooth 
refers  to  the  following  conditions: 

Side  roughness,  1/%-in.  lucite  strips  on  sidewalls  only, 
spaced  on  2-in.  centers  (n  =  0.020) 

Bottom  roughness,  l/4-in.  lucite  strips  on  bottom  only, 
spaced  on  2-in.  centers  (n  =  0.018) 

Smooth,  no  roughness  strips  (n  =  0.011) 

The  indicated  values  of  Manning's  n  were  based  on  measurements  with 
fresh  water  under  steady-flow  conditions.  For  the  series  III  tests. 


Fig.  7*  Comparison  of  experimental 
and  calculated  maximum  velocities 
along  channel 
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table  3  shows  that  the  value  of  m  changes  in  direct  proportion  to  the 
change  in  n  . 

57-  It  is  concluded  that  the  values  of  p  and  k  as  determined 
herein  accurately  reflect  the  effect  of  the  various  independent  quantities 
on  the  tidal  dynamics.  The  length  and  celerity  of  the  tidal  wave  are  also 
given  in  table  3*  Both  X  and  c  show  a  reduction,  in  comparison  with 


the  frictionless  values  X  and  c  , 

o  o 


DISTANCE  ALONG  FLUME  IN  FT 

Fig.  8.  Variation  of  accumulated  rate 
of  energy  dissipation  along  channel 


of  approximately  75  per  cent. 

58.  The  value  of  the  rate 
of  energy  dissipation  (per  unit 
mass  of  fluid),  G  ,  due  to  tidal 
friction  is  calculated  from  equa¬ 
tion  26  and  is  also  given  in 
table  3.  Fig.  8  shows  a  plot  of 
the  rate  of  energy  dissipation  in 
successive  portions  of  the  channel 
divided  by  the  total  rate.  Within 
the  region  of  salinity  intrusion 
this  ratio  is  reasonably  constant. 


Determination  of  Apparent  Diffusion  Coefficient 
and  Boundary  Conditions  at  Ocean  Entrance 


59-  The  fictitious  length  B  marks  the  seaward  extension  of  the 
channel  to  a  point  at  which  the  salinity  is  constant  at  all  times  in  the 
tidal  cycle.  This  distance  is  evaluated  by  assuming  that  the  salinity- 
convection  function  given  by  equation  tO  also  holds  in  this  portion.  It 
is  observed  that  at  the  ocean  entrance  to  the  channel  there  is  a  fairly 
well-defined  time  interval  3  between  the  time  of  low  water  in  the  ocean 
and  the  time  at  which  the  ocean  salinity  first  appears  at  the  entrance  on 
the  flood  tide.  This  is  shown  graphically  in  fig.  9  (test  15)  in  which 
the  variation  of  mean  salinity  with  time  is  shown  for  the  test  station 
nearest  the  ocean  entrance  (x  =  5  ft).  With  the  time  interval  (3  thus 
known,  the  distance  B  is  calculated  from  equation  bO  by  noting  that  the 
salinity  at  low  tide  (t/T  =  0)  at  x  =  -B  is  equal  to  the  salinity  at 
t/T  =  3  at  x  =  +5  ft.  Hence, 
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Fig.  9*  Experimental  salinity  versus  time  curve  for  determination 
of  time  interval  f3  and  minimum  salinity 


(1  +  5B)  e"a/h  =  (1  -  56)  e"a/h(cOS  2nf3) 

from  which, 

B  =  i  (!  -  56)  ea/h(l  -  C0S  2nP) 


(52) 

(53) 


The  constant  W  in  the  low-tide  salinity  distribution  equation  45  can  be 
calculated  by  observing  the  minimum  salinity  at  the  station  nearest  the 
ocean  entrance  (x  =  5  ft),  also  shown  in  fig.  9-  From  equation  4-5, 

(r-)  =  -w(xl  +  B)2 

'  o  min 


hence,  with  x  =  5> 


(5  +  B)2 


and 
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60.  All  of  the  quantities  required  to  calculate  the  instantaneous 
salinity  distributions  given  by  equation  k'J  have  been  determined.  In  addi¬ 
tion,  the  maximum  and  minimum  intrusion  lengths,  equations  5°  and  51,  may 

be  calculated.  For  each  test  the  values  of  0  ,  B  ,  ( s/s  )  .  ,  and  D' 

o  rain  o 

are  tabulated  in  table  4. 

Comparison  of  Experimental  and  Analytical  Salinity  Distributions 

61.  A  graphical  comparison  of  the  calculated  and  experimental  in¬ 
stantaneous  salinity  distributions  at  10  per  cent  time  intervals  throughout 
the  tidal  cycle  is  given  for  test  15  in  fig.  10.  The  agreement,  especially 
for  the  maximum  and  minimum  intrusion  lengths  is  considered  satisfactory. 


Fig.  10.  Comparison  of  experimental  and  analytical  instantaneous  salinity 

distributions  for  a  tidal  cycle 

62.  In  the  following  sequence  of  figures  (figs.  11-15)  the  experi¬ 
mental  and  analytical  salinity  distribution  comparisons  for  other  tests 
are  given  for  both  high-water  (HW)  and  low-water  (LW)  conditions.  In 
these  plots,  the  analytical  salinity  distribution  is  defined  by  the  curves 
and  the  experimental  salinity  distribution  by  the  plotted  points.  Fig.  11 
shows  the  maximum  and  minimum  salinity  distributions  for  tests  l6,  2,  and 
11  (series  i).  This  group  of  tests  shows  the  influence  of  the  fresh-water 


Fig.  11.  Comparison  of  experimental  and  analytical  salinity 
distributions  for  series  I,  tests  l6,  2,  and  11 
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discharge  which  increases  in  the  three  tests  in  proportion  to  1: 2:2.8, 
all  other  conditions  being  constant.  Results  of  the  remaining  tests 
of  series  I  (l8  and  26)  are  shown  in  fig.  12.  Aside  from  a  slight 
roughness  change,  the  primary  variable  in  these  two  tests  is  the  ocean 
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Fig.  12.  Comparison  of  experimental  and  analytical 
salinity  distributions  for  series  I, 
tests  18  and  26 
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salinity  which  is  approximately  five  times  larger  in  test  26  than  in 
test  18. 

63.  The  two  tests  ( 15  and  10)  in  series  II  are  shown  in  fig.  13 
for  different  fresh-water  flow  rates  and  tidal  amplitudes  50  per  cent 


*  IN  FT 


Fig.  13.  Comparison  of  experimental  and  analytical 
salinity  distributions  for  series  II, 
tests  15  and  10 


greater  than  those  of  series  I.  Fig.  l4  gives  a  sequence  of  three  tests 
(l4,  4,  and  9)  in  series  III  (tidal  amplitude  twice  as  large  as  those  of 
series  i)  for  increasing  fresh-water  flow  rates.  Finally,  tests  19  and 
34  in  fig.  15  show  effects  of  ocean  salinity  and  roughness  change  to  com¬ 
plete  the  series  III  tests.  In  all  cases  the  calculated  and  experimental 
salinity  distributions  are  in  reasonable  agreement. 


Fig.  15.  Comparison  of  experimental  and  analytical 
salinity  distributions  for  series  III, 
tests  19  and  34 
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PART  IV:  CORRELATION  OP  DIFFUSION  PARAMETERS 

64.  The  determination  of  the  instantaneous  salinity  distribution  is 

dependent  upon  two  basic  parameters:  the  apparent  diffusion  coefficient 

at  the  channel  entrance  D1  and  the  seaward  excursion  distance  B  .  Both 

o 

quantities  can  be  determined  from  an  experimental  record  of  the  salinity 
variation  over  a  tidal  cycle  at  the  flume  entrance.  However,  for  a  more 
complete  understanding,  a  correlation  between  ,  B  ,  and  the  independent 
quantities  defined  by  equation  1  is  required. 

65.  One  relation  is  immediately  obtained  from  the  Kolmogoroff 

theorem  relating  the  eddy  diffusivity  to  the  rate  of  energy  dissi¬ 

pation  per  unit  mass  of  fluid  G  and  the  mean  eddy  size  £  . 

Dx  =  (constant)  G1/3  (55) 

If  the  scale  of  turbulence  £  is  a  function  of  the  physical  dimensions  of 

the  system  and  is  assumed  to  be  constant,  the  eddy  diffusivity  should  vary 
l/3 

in  proportion  to  G 

66.  Tests  31,  32,  and  33  in  the  WES  flume  were  conducted  using  dyed 
fresh  water  in  place  of  the  saline  ocean  water.  The  dye  concentration  was 
of  the  order  of  1  ppm  and  density  differences  were  therefore  negligible. 

The  coefficients  obtained  are  thus  true  eddy  diffusivities  unaffected  by 
gravitational  convection.  From  table  4  it  is  seen  that  for  tests  32  and 
33  the  ratio  Dq/g1^3  is  equal  to  a  constant  2.8  (Dq  for  test  31  could 
not  be  obtained  due  to  a  lack  of  definition  of  the  experimental  concentra¬ 
tion  curve).  It  is  also  noted  that  the  ratio  D^/g^3  for  all  other  tests 
employing  saline  water  is  greater  than  2.8.  This  ratio  reaches  a  maximum 
value  of  25  for  test  11  which  represents  the  maximum  fresh-water  discharge 
and  maximum  ocean  salinity  for  the  minimum  tidal  range  (i.e.  the  most 
highly  stratified  condition) . 

67.  The  above  observations  are  illustrated  in  fig.  16  which  is  a 

plot  of  and  Dq  versus  the  cube  root  of  the  tidal  energy  dissipation 

rate  G  .  The  zero-density  runs  plot  as  a  straight  line  indicating  that 

l/3 

Dq  is  directly  proportional  to  G  .  In  other  words  the  transfer  of 
salt  upstream  by  turbulent  diffusion  increases  with  the  rate  of  tidal 
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TIDAL  ENERGY  DISSIPATION  RATE,  g'/3 


Fig.  16.  Diffusion  coefficient  versus  rate  of  tidal  energy  dissipation 


energy  dissipation.  The  saline  tests  exhibit  the  opposite  trend,  a  non- 

.  l/3 

linear,  inverse  relation  between  and  G  .  As  the  energy  dissipa¬ 

tion  (and  therefore  the  turbulence)  increases,  the  transfer  of  salt  up¬ 
stream  by  the  combined  processes  of  gravitational  convection  and  diffusion 
decreases.  This  trend  is  indicated  by  the  shaded  region  in  fig.  16  and 
shows  that  as  the  turbulent  mixing  increases,  the  gravitational  convection 
of  salt  is  reduced.  A  measure  of  the  partially  mixed  estuary  is  the  ratio 
of  the  apparent  diffusion  (due  to  the  combined  salt-transfer  proc¬ 

esses)  to  the  turbulent  diffusion  coefficient  Dq  .  Hence,  from  equa¬ 
tion  55 

D'  D' 

o  o 

“o'?75 

1/ 

66.  The  ratio  D^/GA/  is  seen  to  be  a  significant  expression  of 
the  increase  in  salt  transfer  due  to  gravitational  convection.  Hence,  it 
would  be  expected  to  correlate  with  some  function  involving  the  density 
difference  between  the  fresh  water  and  the  ocean.  In  the  MIT  tests  (refer¬ 
ence  10)  such  a  correlation  was  found  between  D^/G'^  and  a  dimensionless 
parameter  of  the  following  form: 


ko 


G  _  rate  of  energy  dissipation  per  unit  mass  of  fluid 
J  “  rate  of  potential  energy  gain  per  unit  mass  of  fluid 


The  gain  in  potential  energy  (per  unit  volume),  as  water  moves  down  the 
estuary  toward  the  sea,  is  due  to  its  increasing  specific  weight,  i.e. 


Ayh 


lary  towari 
( ft -lb  A 
Vft3  J> 


where  Ay  =  difference  in  specific  weight  of  fresh  water  and 


ocean;  thus,  the  gain  in  potential  energy  per  unit  mass  is  therefore 


ATh 

P 


Since  the  average  time  of  retention  of  the  fresh  water  in  the  estuary  is 
tr  =  jj  ,  the  time  average  rate  of  gain  of  potential  energy  per  unit  mass 
of  fluid  is  given  by 


J  = 


(56) 


Hence,  from  equation  26 


G  _  1  c  ( af  M 
J  =  2  u  Vh/  Kyfy 


(57) 


where  the  tidal  wave  celerity  c  =  o/k  and 


M 


_ sinh  2uL _ 

cos  2kL  +  cosh  2pL 


(58) 


The  -values  of  G/j  are  tabulated  in  table  k .  Fig.  17  is  a  plot  of 
D'/G1^  versus  g/j  ,  and  it  is  seen  that  a  significant  correlation  of  the 
apparent  diffusion  coefficient  has  been  obtained  in  terms  of  the  tidal  and 
fresh-water  flow  characteristics  of  the  estuary.  The  parameter  G/j  is 
analogous  to  a  Richardson  number  which  expresses  the  stability  criterion 
for  flows  with  vertical  density  gradients.  The  parameter  is  furthermore  a 
quantitative  measure  of  the  degree  of  stratification  or  mixing  in  an 
estuary.  Increasing  values  of  G/j  Indicate  increasingly  well -mixed  con¬ 
ditions;  for  the  case  of  zero  gravitational  effects,  G/j  approaches 
infinity. 

69.  The  remaining  problem  in  the  one -dimensional  treatment  is  to 
find  a  correlation  between  the  seaward  excursion  (defined  by  the  fictitious 
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Fig.  17.  Correlation  of  apparent  diffusion  coefficient  D' 
with  stratification  number  ° 


length  B)  and  the  independent  quantities  of  the  system.  Fig.  9  shows  the 
time  variation  of  the  salinity  at  the  ocean  entrance  of  the  flume.  It  is 
noted  that  full  ocean  salinity  exists  at  this  point  during  approximately 
one-half  the  tidal  cycle  and  that  during  the  other  half  of  the  cycle  the 
average  salinity  is  reduced.  The  physical  explanation  is  that  (during  the 
major  portion  of  the  flood  tide)  the  fresh-water  inflow  is  stored  in  the 
estuary.  On  the  ebb  cycle  the  removal  of  the  fresh  water  causes  a  dilu¬ 
tion  as  indicated  by  the  observed  salinity  curve.  It  is  apparent  from 
equation  53  that  the  time  lag  P  and  therefore  the  parameter  B  are  de¬ 
pendent  on  the  shape  of  the  salinity  versus  time  curve  at  the  flume 
entrance. 

70.  The  shape  of  this  curve  is  given  by  equation  kj;  with  x  =  0  , 


si0**) 

8 


■f-  - - - 5  [l  -  ea/h<1  -  °°a  .  8b12  \ 

l  2D'  B5  L  J  J 
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If  a/h  is  small  compared  to  unity,  the  exponential  term 


a/h( 1  -  cos  it)  ,  a  ,  . 

e  '  sr  1  »-  -  ( 1  -  cos  at) 

and  since  8  =  ,  equation  59  may  be  written 

o 


=  exp 


UB  |  ,  Uo  n  v  2 

2D"5"  I  1  ”  oB  1  C°S  rt)  J 

o  ' 


The  salinity  variation  vith  time  at  the  estuary  entrance  is  therefore  seen 

UB  u0 

to  be  a  function  of  two  dimensionless  parameters  — r  and  ~  .  The  phys- 

Dq  cri3 

ical  significance  of  the  latter  quantity  may  be  given  as  follows:  The 


tidal  excursion  is  the  average  longitudinal  distance  traveled  by  a  parti- 

u0  T 

cle  in  a  tidal  cycle;  this  distance  is  proportional  to  — r —  and  since 

/  uo  2 

a  =  2tr/T  ,  tidal  excursion - . 

u0  0 

71 •  Therefore,  —  is  a  parameter  proportional  to  the  ratio  of  the 
tidal  excursion  to  the  seaward  excursion  of  salinity  B  .  It  has  been 
found  that  this  quantity  is  related  to  the  stratification  parameter  G/j 


previously  discussed. 

72.  Fig.  18  is  a  plot  of  uB/uq  versus  G/j  for  all  of  the  salinity 


STRATIFICATION  NUMBER,^ 

Fig.  18.  Correlation  of  seaward  excursion  distance  B 
with  stratification  number 


tests  tabulated  in  table  4.  It  is  concluded  that  a  reasonable  correlation 
exists  and  that  well-mixed  conditions  are  characterized  by  small  values  of 
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PART  V:  CONCLUSIONS 

Longitudinal  Salinity  Distribution  in  a  Uniform  Estuary 

73*  It  is  concluded  that  the  objectives  of  the  one -dimensional  phase 
of  this  investigation  have  been  accomplished.  Analytical  expressions  have 
been  developed  to  describe  the  instantaneous  distribution  of  the  mean  ver¬ 
tical  salinity  throughout  the  tidal  cycle,  and  thus  the  maximum  and  minimum 
salinity-intrusion  lengths  for  the  uniform  estuary  are  also  specified. 

74.  In  general,  the  following  independent  quantities  must  be  speci¬ 
fied  initially:  estuary  length,  mean  depth,  ocean  tidal  amplitude,  period, 
ocean  salinity  concentration,  and  fresh-water  discharge.  The  salinity  dis¬ 
tribution  equations  involve  five  basic  quantities  in  addition  to  those 
specified  above:  the  tidal  damping  modulus  q  ,  the  tidal  wave  number 
k  =  2rt/l  ,  the  maximum  tidal  velocity  at  the  ocean  entrance  uq  ,  the  ap¬ 
parent  diffusion  coefficient  at  the  ocean  entrance  ,  and  the  seaward 
excursion  B  .  The  first  three  (p  ,  k  ,  uq),  being  exclusively  tidal  pa¬ 
rameters,  are  determined  from  observations  of  the  tidal  elevations  within 
the  estuary  by  means  of  the  damped  co-oscillating  tidal  theory  developed 
in  paragraphs  18-32.  This  method  of  determining  the  tidal  characteristics 
from  elevations  (which  in  the  case  of  actual  estuaries  are  either  already 
known  or  readily  measured)  shows  promise  for  future  field  studies. 

75*  It  has  been  demonstrated  that  the  remaining  quantities  and 
B  can  be  determined  from  an  observation  of  the  salinity  variation  with 
time  at  the  ocean  entrance  as  shown  in  paragraphs  59-60.  On  the  basis  of 
these  determinations  significant  correlations  have  been  developed  between 
,  B  ,  and  the  tidal  and  physical  specifications  of  the  estuary.  In  par¬ 
ticular,  and  B  are  correlated  with  G/j  ,  the  ratio  of  the  mean  rate 

of  turbulent  energy  dissipation  to  the  rate  of  increase  of  potential  energy 
(due  to  the  change  in  salinity  along  the  estuary).  It  is  suggested  that 
this  ratio  be  called  the  stratification  number  since  it  has  been  shown  to 
be  a  quantitative,  measure  of  the  extent  of  vertical  mixing  in  an  estuary. 

76.  If  use  is  made  of  the  approximation  e-2a/k  =  1  -  2a/h  ,  the 
equation  (49)  for  the  maximum  salinity-intrusion  length  may  be  written: 


2a 

h 


L 


t 
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(49a) 


From  the  above  it  is  seen  that  the  intrusion  length  in  a  uniform  estuary 
varies  in  a  relatively  complex  manner  with  the  various  tidal  and  diffusion 
quantities.  For  the  flume  tests  a  quantitative  evaluation  of  the  influ¬ 
ence  of  these  quantities  on  the  intrusion  length  may  be  made  by  comparing 

the  values  of  L,  indicated  in  figs.  11-15- 
x> 

77.  Equation  49a  may  also  be  written  in  dimensionless  form  by  using 
the  tidal  excursion  parameter  u  /o  as  a  reference  length;  thus. 


The  ratio  of  the  maximum  salinity-intrusion  distance  to  the  tidal  excur¬ 
sion  therefore  depends  on  three  dimensionless  quantities: 

a.  a/h  ,  the  ratio  of  ocean  tidal  amplitude  to  estuary  mean 

depth 

b.  aB/u  ,  the  ratio  of  the  seaward  salinity  excursion  to  the 

tidal  excursion 

c.  B'/UB  ,  the  ratio  of  the  diffusion  distance  D'/lJ  to  the 

seaward  salinity  excursion 

The  quantities  and  B  are  given  in  terms  of  the  stratification 

number  G/j  in  figs.  17  and  l8.  The  stratification  number  can  be  evalu¬ 
ated  from  the  tidal  parameters  a,  0  ,  k  ,  and  p  ;  the  estuary  length 
L  and  mean  depth  h  ;  the  mean  fresh-water  flow  rate  U  ;  and  the  ocean 
water-fresh  water  specific  weight  difference  A7/7  as  given  by  equation  57* 
78.  On  the  basis  of  the  specific  information  presented  in  this 
study,  it  is  possible  to  make  quantitative  predictions  regarding  salinity 
intrusions  in  partially  or  well-mixed  estuaries  of  essentially  uniform 
cross  section  without  prior  knowledge  of  any  existing  salinity  conditions. 
In  addition,  the  factors  which  cause  changes  in  salinity  intrusions  have, 
for  the  first  time,  been  quantitatively  evaluated.  Thus  by  means  of  the 
intrusion  equations,  the  effect  of  changes  in  the  fresh-water  discharge, 
channel  depth,  etc.,  can  be  predicted.  Of  even  greater  importance  is  the 
possibility  that  the  diffusion  and  tidal  parameters  developed  as  a  result 
of  this  study  will  have  important  significance  in  understanding  the  in¬ 
trusion  mechanics  of  estuaries  of  nonuniform  geometry. 
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Significance  of  the  Stratification  Number 


79-  Previous  investigators  of  the  problem  of  salinity  intrusion  in 
tidal  estuaries  have  proposed  that  the  ratio  of  the  volume  of  fresh-water 
inflow  in  a  tidal  cycle  to  the  tidal  prism  is  an  important  parameter  de¬ 
scribing  the  mixing  characteristics  of  the  estuary.  The  tidal  prism  is  de¬ 
fined  differently  by  various  writers.  In  some  cases  it  is  defined  as  the 
change  in  the  volume  of  the  estuary  between  low  water  and  high  water.  In 
other  cases  it  has  been  taken  equal  to  the  volume  of  sea  water  entering 
the  estuary  on  the  average  flood  tide.  The  latter  definition  differs  from 
the  former  in  that  it  does  not  consider  the  volume  of  fresh-water  inflow 
during  the  flood  cycle  as  part  of  the  tidal  prism.  Another  ambiguity 
arises  as  to  the  meaning  of  the  tidal  prism  where  (as  is  almost  always  the 
case)  the  time  of  high  water  or  low  water  occurs  at  different  times  along 
the  estuary.  According  to  the  first  definition,  the  tidal  prism  would  be 
the  volume  contained  in  the  envelope  between  high-  and  low -water  levels. 

The  objective  of  this  discussion  is  to  relate  the  stratification  number  to 
the  fresh  water-tidal  prism  ratio.  For  this  purpose  the  tidal  prism  is 
taken  as  the  volume  of  sea  water  entering  the  estuary  on  the  flood  tide. 

If  the  tidal  velocity  at  the  estuary  entrance  varies  harmonically  with 
time,  the  average  velocity  over  a  half  cycle  (flood  tide)  is  2/m  times 
the  maximum;  hence,  the  tidal  prism  P  is 


2  T  un  bhT 

P .  =  -  u  (bh)  |  =  -2 - 

t  «  o  '  '  2  n 


The  volume  of  fresh-water  inflow  during  a  full  tidal  cycle  is 


Qf  T  =  UbhT 

Hence,  the  ratio 

V;  su 

pt  =u0 

This  ratio  is  plotted  against  the  stratification  number  in  fig.  I9¬ 
60.  From  equation  57  it  can  be  seen  that  for  a  given  estuary  in 
which  the  only  variable  is  the  fresh-water  velocity  U  ,  the  stratification 
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Fig.  19-  Relation  between  fresh-water  volume: tidal  prism  ratio 
and  stratification  number 


number  G/j  is  proportional  to  l/u  .  The  tidal  prism  ratio  under  these 
conditions  is  proportional  to  U  ;  therefore,  it  would  be  expected  that 
for  a  sequence  of  runs  in  which  only  U  is  varied,  the  relation  between 
the  two  parameters  would  be  linear  in  logarithmic  coordinates.  Fig.  19 
shows  the  expected  result  for  the  test  sequences  (11,  2,  16),  (10,  15), 
and  (9,  14).  The  important  conclusion  is  that  the  tidal  prism  ratio  is 

not  indicative  of  the  extent  of  mixing  if  any  other  quantity  such  as  tidal 
amplitude,  depth,  or  roughness  is  a  variable.  In  other  words,  equal  values 
of  the  tidal  prism  ratio  do  not  necessarily  indicate  similar  mixing  condi¬ 
tions;  however,  equal  values  of  the  stratification  number  do  indicate 
similar  conditions.  This  conclusion  is  demonstrated  by  fig.  20  which 
shows  on  a  dimensionless  basis  the  vertical  salinity  gradients  for  tests 
11,  10,  1 6,  and  l4.  The  plotted  profiles  are  for  a  station  40  ft  from  the 
ocean  tidal  basin  and  correspond  to  the  conditions  existing  at  approxi¬ 
mately  one-quarter  of  a  tidal  cycle  after  the  time  of  high  water  in  the 
basin,  at  which  time  the  vertical  salinity  gradients  for  all  runs  are  a 
maximum.  Fig.  20  clearly  indicates  the  correlation  between  the  stratifica¬ 
tion  number  and  the  vertical  salinity  gradient.  Test  11  (G/j  =  20)  is  the 
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Fig.  20.  Vertical  salinity  gradients  in  relation 
to  stratification  number 

most  stratified  and  test  l1)-  (G/j  =  202)  the  least  stratified  of  all  the 
runs  with  30  ppt  salt  in  the  ocean  basin.  Of  particular  interest  here  is 
that  tests  10  and  16  are  from  different  series,  with  test  1 6  having  two- 
thirds  the  tidal  amplitude  and  one-third  the  fresh-water  flow  rate  of  test 
10.  The  stratification  numbers  for  these  two  tests  are  approximately  the 
same,  and  the  vertical  salinity  gradients  as  shown  in  fig.  20  are  also 
similar.  The  values  of  the  tidal  prism  ratio,  however,  differ  by  a  factor 
of  two  as  is  readily  verified  from  fig.  19 • 

8l.  It  is  concluded  that  the  stratification  number  is  sensitive  to 
the  various  factors  affecting  the  salinity  regime  in  a  uniform  estuary. 

The  tidal  prism  ratio  is  sensitive  only  to  changes  in  fresh-water  flow, 
all  other  conditions  being  equal. 


PART  VI:  RECOMMENDATIONS 


82.  It  is  proposed  that  the  general  investigation  of  salinity  in¬ 
trusion  be  continued  along  the  lines  of  the  four -phase  program  originally 
outlined  by  the  Tidal  Hydraulics  Committee  (see  Preface).  It  is  suggested 
that  the  second  and  third  phases,  i.e.  the  vertical  mixing  of  fresh  and 
salt  water,  the  resulting  vertical  salinity  distribution,  and  the  vertical 
distribution  of  current  velocities  as  affected  by  salinity  distribution,  be 
treated  as  a  single  unit.  In  essence,  the  second  and  third  phases  imply  a 
two-dimensional  analysis  of  the  salinity-intrusion  problem,  and  they  are  a 
logical  extension  of  the  one -dimensional  analysis  treated  in  this  report. 

83.  Recent  research  in  the  MIT  Hydrodynamics  Laboratory  under  non- 
tidal  conditions  (reference  7)  and  earlier  field  investigations  by  Prit¬ 
chard  in  the  James  River  have  shown  the  feasibility  of  a  two-dimensional 
analysis.  Vertical  velocity  distributions  can  be  derived  from  horizontal 
and  vertical  salinity  gradients  through  the  conservation-of-mass  equation. 
The  time  average  vertical  velocity  distributions  are  of  the  familiar  "s" 
shape  with  a  seaward  flow  near  the  surface  and  a  landward  flow  near  the 
bottom.  The  MIT  tests,  for  a  range  of  partially  mixed  conditions,  showed 
surprisingly  large  magnitudes  of  the  ratio  of  landward  bottom  velocity  to 
average  vertical  velocity.  Of  particular  interest  is  the  fact  that  a  rela¬ 
tion  seems  to  exist  between  this  local  bottom  velocity  ratio  and  the  local 
densimetric  Froude  number.  The  densimetric  Froude  number  is  based  on  the 
fresh-water  velocity,  depth,  and  local  density  difference  ^x/p  •  In  the 
landward  direction,  this  quantity  increases  (since  Apx/p  is  in  the  de¬ 
nominator)  and  the  bottom  velocity  decreases.  These  results  are  obtained 
in  the  case  of  a  Bteady-state  salinity  distribution  in  homogeneous  turbu¬ 
lence.  The  experimental  data  already  obtained  in  the  WES  flume  can  be 
analyzed  in  a  similar  manner.  More  specific  information  on  variations  of 
the  landward  bottom  velocity  or  flow  predominance  along  an  estuary  should 
be  useful  in  view  of  its  demonstrated  importance  to  sediment  movement  and 
deposition.  This  leads  directly  to  the  fourth  phase  of  the  general  study. 

84.  As  a  parallel  investigation  it  is  also  recommended  that  the  one¬ 
dimensional  approach,  as  contained  in  this  report,  be  broadened  to  include 
a  study  of  the  longitudinal  distribution  of  salinity  in  an  estuary  of 
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nonuniform  geometry.  For  this  purpose  available  test  data  from  the  modi 
fied  Delaware  estuary  model  could  be  utilized  to  advantage. 
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Table  1 


Summary  of  Flume  Tests  Used  in  Analysis 

Flume  depth  at  msl  =  0.5  ft  Tidal  period  =  600  sec  (except  30- 3) 

Length  of  flume  =  327  ft  Flume  width  =  0.75  ft 
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30 
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Side 
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Tidal  Analysis,  Test  Sj 
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Table  3 

of  Tidal  Characteristics 


General  Summary 

Mean  depth,  h  =  0-5  ft  Width,  b  =  0.75  ft 

Tidal  period,  T  =  600  sec  (except  where  noted)  cq  =  -J gh  =  1.0  ft/ sec 

0  =  0.010  l/sec  X  =  c  T  =  2l00  ft 

'  A  ft 


Length,  L  =  327  ft 

k  = 
0 

0.0026 

radian/ ft 

Test 

8o 

U 

k 

u 

X  -  SI 

O 
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h3|>" 

G  X  103 

U 

0 
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No. 

Roughness 

PPt 
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radian/ft  radian/ft 

ft 
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l/ft 

Series  I, 

Tidal  Amplitude  =  0, 

.05  ft 

29 

Side 

0 

0 
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0.0011 
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3.18 

— 

0.11 

— 

30-3* 

Side 

0 

0 
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0.0021 
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-- 

0.18 

-- 
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Side 
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0.02 
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15 

Side 

29.5 

0.02 

0.0035 

0.0017 

1790 

3.00 

2.20 

0-59 

0.0026 

10 

Side 

26.8 

O.O56 

(T15) 

(T15) 

(1790) 

(3-00) 

(2.20) 

(0.59) 

(0.0026) 

32 

Side 

Dye 

0.02 

(T15) 

(T15) 

(1790) 

(3-00) 

(2.20) 

(0-59) 

(0.0026) 

Series  III, 

,  Tidal  Amplitude  =  0.10  ft 
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— 
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of -Diffusion  Parameters  for  Salinity  Intrusion 
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U.  S.  Army,  Committee  on  Tidal  Hydraulics,  CE,  ONE-DIMENSIONAL 
ANALYSIS  OF  SALINITY  INTRUSION  IN  ESTUARIES,  by  A.  T.  Ippen  and 
D.  R.  F.  Harleman.  June  1961,  x,  52  pp  -  illus  -  tables.  (Tech¬ 
nical  Bulletin  No.  5) 

Unclassified  report 

An  experimental  and  analytical  study  of  the  basic  factors  which  de¬ 
termine  the  instantaneous  longitudinal  distribution  of  salinity  in  a 
partially  mixed  estuary  of  constant  section  is  presented.  The  mass- 
transfer  equation  is  solved  in  two  parts:  (a)  the  quasi-steady- 
state  convective-diffusion  problem  for  an  observer  moving  with  the 
tidal  velocity,  and  (b)  the  cyclic  translation  of  saline  water  due 
to  the  tidal  velocity.  In  the  steady-state  portion  the  apparent 
longitudinal  diffusion  coefficient  is  a  function  of  the  turbulence 
induced  by  the  tide  and  of  the  internal  circulations  Induced  by  the 
density  difference.  An  expression  is  developed  for  the  intrusion 
length  as  a  function  of  the  estuary  length,  mean  depth,  tidal  ampli¬ 
tude  and  period,  fresh-water  discharge,  ocean  salinity,  and  estuary 
roughness.  A  stratification  number  is  defined  which  expresses  the 
degree  of  stratification  or  mixing  in  an  estuary. 
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